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Abstract 

We study the existence and properties of the non-equilibrium steady state 
which arises by putting two copies of systems at different temperatures into 
a thermal contact. We solve the problem for the relativistic systems that are 
described by the energy-momentum of a perfect hydro with general equation 
of state (EOS). In particular, we examine several simple examples: a hydro 
with a linear EOS, a holographic GET perturbed by a relevant operator and 
a barotropic fluid, i.e., P = P{£). Our studies suggest that the formation of 
steady state is a universal result of the hydro regime regardless of the kind of 
fluid. 
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1 Introduction 


Although the study of equilibrium states has been the focus of most research in 
many body systems and condensed matter physics, much of the interesting phenom¬ 
ena around us are far from equilibrium. However, thermodynamic study of non¬ 
equilibrium states are less advanced since the dynamical equations governing ther¬ 
modynamic variables are highly non-linear differential equations. Unlike the thermo¬ 
dynamic equilibrium where the processes are reversible and time independent, for a 
system in a non-equilibrium state processes are generally irreversible and time de¬ 
pendant. That’s what makes the study of systems far from equilibrium so difficult. 
Nevertheless, waiting long enough, most non-equilibrium systems tend to approach 
to the state of thermodynamic equilibrium unless there is a continuous flow of matter 
or energy to or from the system. 

Yet an interesting subclass of non-equilibrium systems is the steady state which 
is reached by a system subject to a constant driving force. In a non-equilibrium 
steady state (NESS) there is no time variation, i.e., all thermodynamic properties 
are time independent, although the entropy production is non-zero and the system 
is recognised by the presence of fluxes. The fluxes are the flows of currents of the 
conserved variables driven by the gradient of their corresponding thermodynamic 
conjugate, e.g., the energy current is driven by the temperature gradient, and so 
on. Since the dynamics of such a system is dissipative, in order to maintain the 
steady state one needs the injection of energy at one boundary and subtraction at 
the other. One simple example of such NESS is a system confined between two heat 
baths (reservoirs) at different but constant temperature, transporting heat from one 
to another by the energy current which is generally proportional to the temperature 
gradient. 

However in a recent studies of one dimensional quantum critical systems described 
by conformal field theory (GET), Bernard and Doyon pQ found an interesting steady 
state with a nontrivial energy current choosing a different framework; Instead of 
putting a system in contact with external reservoirs, they picked two copies of some 
large quantum systems at different thermal equilibrium, e.g., different temperatures, 
and then glued them together at a contact point which result in the energy transfer 
from one to another. After long enough time, a steady state would be established 
and parts of each system far from the contact point would effectively behave like heat 
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baths. Another interesting point abont their set-np is that althongh the temperatnre 
prohle is flat, there exists a non-zero energy cnrrent. That is, nnlike the nsnal case 
for steady states, the energy cnrrent is not simply a fnnction of temperatnre gradient. 
Moreover, the nniversal character of the resnlting steady state is noticeable; the energy 
transfer only depends on the nniversal constants and the temperatnres of the initial 
copies they start with. 

Motivated by these resnlts, the existence of snch a steady state and its nniversality 
has been investigated for CFT in higher dimensions [2lE]. Therefore, one legitimate 
qnestion might be whether this is a particular feature of CFT, or one may get similar 
answer starting from initial systems with more general equation of state. In this 
manuscript we are trying to explore this question relying on the holographic insight 
given in [3]|i] there it has been shown that the non-equilibrium steady state on the 
conformal boundary is in fact dual to the Lorentz boosted black brane in the bulk. 
According to the AdS/CFT correspondence, the stress tensor of the steady state could 
be obtained from the metric of the dual gravity theory in the bulk. Therefore, the 
steady state is completely specified if one obtains its temperature and boost velocity, 
i.e., {T, n}, in terms of the temperatures of the initial systems at the instant of 
contact, denoted by Ti and Tj. for one system on the left and the other on the right. 
In fact, at sufficiently large scales two initial systems look like asymptotic reservoirs 
with the steady state as the intermediate state interpolating between two reservoirs 
on the left and right. Then in this framework, the problem is similar to a Riemann 
problem with initial boundary condition Ti for x < 0 and for x > 0, where a 
possible solution would be two shock waves emitting at the point of contact and 
moving to the opposite directions. Note that the development of the steady state in 
this set-up also requires that two shock waves propagate in the opposite directions 
without splitting and decomposition, for if this happens the composite waves could 
lead to cascades and finally thermalizing the system. 

Riemann problem for classical and relativistic hydrodynamics with various equa¬ 
tion of states is widely studied in the literature, see for example HI El El El i. We 
follow similar approach of solving Riemann problem to study NESS while considering 

^We prefer to follow the method introduced in [3] as opposed to the ansatz given in [5]. The 
former, due to holographic analysis, has less parameters and is easier to deal with particularly if one 
interested in studying a more general equation of state like barotropic case that we investigate in 
this paper. 
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to have a relativistic perfect fluids for the left and right hydro before bringing them 
into a thermal contact. We solve the problem assuming a general equation of state, 
e.g., barotropic fluid. In general, we conclude that the formation of steady state after 
thermal contact has nothing to do with conformal held theory or even with integrable 
models; it is a universal result in the hydro regime for just about any kind of huid. 

Therefore, the outline of this paper is as follows: in section H] we briehy review 
shock wave solutions to the Riemann problem as well as Rankine-Hugoniot jump 
condition which are essential in our calculations of specifying the steady state. Then 
in section [3], starting with the perfect huid assumption for the initial systems at 
the instant of contact, we investigate the steady state properties for the systems with 
more general equation of state (EOS) than the conformal huid. In particular, we study 
three examples: systems with linear equation of state, the QFT model described in [9] 
as a small deviation from GET, and hnally systems with barotropic equation of state 
where pressure is only a function of energy density however we restrict our calculation 
to the case where the temperature diherence between two systems are very small. We 
close the paper by some discussion about limitation on the stability and the existence 
of shock solutions and also give a brief conclusion in section HJ 

2 Shock waves and jump condition 

Consider a relativistic huid described by the energy density E and the pressure P, 
and then introduce a small perturbations to the pressure and density of the system. 
It is straightforward to show that these perturbations, the so-called sound waves’^ 
propagate through the huid with the velocity Cg given by 

9 dP 

Is' 

known as the speed of sound' in literature. If we assume that the sound waves have 
inhnitesimal amplitude, i.e., the perturbations are very small, then the speed of sound 
is nearly uniform throughout the huid. However, in general the speed of sound is a 
function of density which result in the crest of wave to move faster than the trough. 
When the crest overtakes the trough a shock wave can form due to the steepening 
of wave. A shock front is a surface that marks a sudden jump in the density and 
pressure of the huid. Although this is a way that shocks could form, in general, 
shock waves are characterized by a rapid, discontinuous change in the density and 
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pressure of the system. The shock is bounded into an inhnitesimal region where the 
fluid properties, such as density and pressure, immediately before and after being 
shocked are linked by the jump conditions. These conditions, usually referred to as 
the Rankine-Hugoniot jump conditions, are derived from the conservation laws. 

Consider a conserved quantity q{x,t) in (l+l)-dimensionJ§ satisfying the hyper¬ 
bolic conservation law HE], he.. 


dtq + d^f{q) = 0, 


( 2 ) 


where / represents the flux of q. Now assume that there exists a solution q which 
has a discontinuity, i.e., a shock, along the curve x = ^(f); this is known as the weak 
solution of the PDE ([2]) [10]. Then choose the interval < x < such that it 
intersects with x = ^{t) at time t. Integrating equation ([2|) over this interval, yields 


'■CW 

If we apply Leibnitz rule to 



ri+ \ /■?+ 

q{x,t)dx + / q{x,t)dx \ + 

Jm 


d^f{q{x,t))dx = 0. 


using the fact that the conservation equation 
satished on either side of the discontinuity, we obtain 

di{t) 


dt 


-{qi - Qr) = fl- fr, 


(3) 


IS 


(4) 


where qi and denote the values of q on the left and right sides of x = ^(f), respec¬ 
tively, and 

fl = f{qi) , fr = f{qr) • 

Equation (j3D is called the Rankine-Hugoniot jump condition and is usually written in 
shorthand notation as 

Us[q] = [/] , (5) 

where [q] = qi — qr and [/] = fi — fr are the jumps across the discontinuity (or shock) 
and Us = df(f) /dt is referred to as the shock speed. 

In more physical systems there are more than one conservation law and the dy¬ 
namical behavior of the fluid is governed by a system of hyperbolic conservation laws 
as 

dtQ + d,F = 0, ( 6 ) 


^The hyperbolic conservation laws originally studied in HIS] for a theory in (l-l-l)-dimensions. 
However, one can follow a similar approach for a theory in higher dimensions by imposing a symmetry 
so that nothing happens in the remaining transverse directions. 
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where Q = (gi,..., and F is a function of gi,..., Therefore, one would derive 
n jump conditions corresponding to each equation in (jH]) as 


ti.lQ] = [F]. 


m 


However, to obtain a full description of the fluid, we also need to take into account the 
equation of state (EOS) which relates different thermodynamic variables of the sys¬ 
tem. Although the conservation laws and the jump conditions are valid for any EOS, 
the deterministic role of EOS on the nature of the shock waves and their propagations 
has been widely considered in literature, see for example mi El [13]. 

Therefore, a key concern in studying the fluid dynamics is the existence and 
uniqueness of the shock solutions in a given fluid, the so-called Riemann problem. In 
general, the Riemann problem is the initial value problem for a system of conserva¬ 
tion laws where the initial data are a pair of constant states separated by a jump 
discontinuity at x = 0, i.e.. 


Qoix) 


Qi if X < 0, 

Qr if X > 0. 


( 8 ) 


In the next section, we study steady state solutions for the perfect fluid with vari¬ 
ous EOSs and examine the constraints to have the solution to the Riemann problem. 


3 Perfect fluid and steady state solution 

Consider two isolated semi-infinite systems each described by a relativistic perfect 
fluid, i.e., 

= (p + ^ (9) 

and governed by the conservation equation 

= 0 . ( 10 ) 

Let’s denote the temperature of one system by T; and the other by and assume 
that the energies of two systems are just a function of temperature, i.e., Si^r = £iTi^r)- 
Now we want to bring two systems into contact at f = 0 and study the possible 
solutions. Indeed, this is a Riemann problem with initial values T; for x < 0 and 
Tr for X > 0. As mentioned in the introduction, this problem has been recently 
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studied for conformal fluid ID Eli- Our goal is to study the problem for an arbitrary 
equation of state while we follow an approach adopted in [3] to examine a possible 
development of the steady state through the propagation of shock waves emanating 
from the contact point. 

It is more evident that the conservation law (fTOj) is a system of hyperbolic equa¬ 
tions when it is written as follows 


dtT^^ + d^T^^ = 0 , 

dtT^^ + d^T^^ = 0 . 

Therefore comparing with ([6]) we get 



(11) 


( 12 ) 


Then all we need to solve for the shock solution are the Rankine-Hugoniot jump 
conditions: 

UsiQl — Qr) = Fi — Fr (13) 

where Q and F are given by flT^ with Qi > Qr and Fi> Fr. 

The existence and uniqueness of the shock solution for an arbitrary EOS is highly 
influenced by the thermodynamic properties of the fluid as well as the nature of the 
conservation laws. Nevertheless, for a relativistic perfect fluid satisfying the hyper¬ 
bolic laws ffTTj) . the existence and uniqueness of the shock solutions has been shown in 
[7] provided the difference between quantities on the left and right side of discontinuity 
is small and the following relations are fulhlled: 


Cg 1, 


<PP ^ (l-cg)cg 
d£^ - P + S 


(14) 


The hrst inequality ensures the physical requirement that the speed of sound is less 
than the speed of light, where we have set c = 1 throughout this paper. 

In the following we study the steady state solutions considering three types of 
equation of state: a) linear EOS, i.e., P = aS, b) the EOS for a {d + l)-dimensional 
QFT introduced in [9] which is in fact a holographic GET perturbed by a relevant 
operator, and c) general barotropic EOS, i.e., P = P{£) where we assume the left 
and right quantity are slightly different, i.e., £i — £r ^ £i,r ■ 
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3.1 Linear EOS 


Consider we start with two sem-infinite fluids described by a linear EOS: 


P = aS, 


( 15 ) 


where from the definition © we infer a = c^. Next we bring two systems into 
instantaneous thermal contact aX x,t = 0. We show that at long enough time the 
steady state forms and we obtain the properties of that. In fact what we have, looks 
like a Riemann problem with initial values 


£o{x) 


El if X < 0 , 

8r if X > 0 . 


(16) 


A consistent solution has two shock waves propagating in opposite direction with 
respect to each other and therefore all we need to solve are the jump conditions fll3p 
for EOS fll5p . At long enough time we have three regions: 

1 . region left described by an ideal fluid with energy momentum tensor: 


rjnflV' 


[ El 0 0 

0 aEi 0 


0 ^ 
0 


y 0 0 0 aEi y 


(17) 


2 . steady state region in the middle described by a boosted fluid with energy E 
and the boost velocity v: 


rr = aEi + (cr-i + , 


(18) 


where = diag(— 1 , 1 , • • • , 1 ) and = 7 ( 1 , n, 0 , • • • , 0 ) with 7 = l/\/l — v^. 

3. region right described by an ideal fluid with energy momentum tensor: 

/Er 0 0 0 \ 

0 aEr 0 0 

rpfll/ _ ' 

r 

(^0 0 0 aEr j 


(19) 


Note that Ei> E > Er- Now we apply the jump relation flT^ for two shock waves: 







I. One shock wave is moving with the speed ui to the negative x direction from 
steady state region with energy £ to the left region with energy £i > £, therefore we 
have: 


uiiQi — Qs) — Fi — Fs, ( 20 ) 

where Q and F are given in flT^ . 

II. Another shock wave is moving with the speed Ur to the positive x direction 
from steady state region with energy £ to the right region with temperature £r < 
therefore we have: 


UriQs - Qr) = Fg- Fr . ( 21 ) 

combining equations (120]) and (1211) immediately gives: 

n~'xx _ n~'xx 

rjntx _ r 

til Xlf 

where we have used flTTD and (ITOD . Moreover, substituting (1171118^ into the equations 
(120112 ip . we obtains the four unknowns in terms of the known boundary 

conditions, i.e., £r{x > 0) and £i{x < 0) as: 


= a 


£l - £r 

til -\- Uf 


( 22 ) 



Ur = 


X + tx 
X + fx-i 
X- 1 


V(x + cx i)(x + a) 


(23) 

(24) 


where x = \f£\j£r- One should note that a = < 1 to fulhl relativistic constraint. 

Accordingly the second expression in (fT4p is satisfied since 0 < a < 1 and therefore 
the above shock solution is the only solution to this set up for all £i and £’,-1^ Moreover, 
for the special case of a = l/d, where d is the spatial dimensions, the equations (1221 
[24p reduced to the one in [5] for conformal fluid considering the fact that the energy 
of the conformal fluid is related to the temperature as £( oc thus x here coincide 
with X in that reference. 

^It is proved that for a linear EOS the uniqueness constraint (ini) is valid for all Ei and £r while 
for general EOS the constraint is valid provided left and right quantities are close. We refer the 
interested reader to [7] for the proof. 
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3.2 A QFT model 


We start with a ((i+l)-dimensional QFT in the regions left and right where the action 
is given by 

Sqft = Soft ^ J d'^^^xO{x) (25) 

where A has dimension d + 1 — A and the dimensionless qnantity -C 1. 

Note that the unitarity bonnd for a scalar operator introduces a lower bound on the 
conformal dimension as A > (d —1)/2, for most recent study see [H]. Also the action 
fl25|l describes a holographic CFT perturbed by a relevant operator which requires 
A < d + 1 and in the following we investigate the development of the steady state for 
such a fluid. Recently in [T5], a similar study has been done for a (l+l)-dimensional 
CFT which is perturbed by an irrelevant operator. 

The energy density and pressure of a fluid described by fl25ll at finite temperature 
has been studied perturbatively in [9] and is given by: 

£(r) = 1-a (26) 

■P(r) = Arm 1^1 _ +..., (27) 


where a = (2A — d — 2)/d and A is proportional to the central charge of CFT. 

Using ([1]) the speed of sound up to second order in perturbative expansion is given 
by 

2(d+1 - 2A)(d+1 - A) 


2 _ 2 
^s,QFT — ^s,CFT 


1 + 


A^ 


( 28 ) 


rf((i+l) p(d+l-A) 

where c^s,cft = ^/d. The second term in the square bracket is negative for (d+l)/2 < 
A < d + 1 indicating that c^s,qft < ^\cft- This is indeed a reasonable expectation: 
in the study of QCD thermodynamics in four dimensions, it is known that in the hot 
QCD plasma the speed of sound is approaching the Stefan-Boltzmann limit, equal to 
conformal limit, i.e., sb — 1/3, from below [I6]. For A < (d + l )/2 it is evident 
that AgQpT > ^^sCFTi violating the Stefan-Boltzman limit, yet an allowed solution 
by unitarity considerations. Also, using the expression fl28p one can compare the 
speed of sound in the left and right fluid. For the conformal dimension in the range 
< A < d + 1 the second term in the bracket is negative, therefore assuming 
Ti > Tr one easily conclude > c^s^Qppp which supports our intuition 
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that the speed of sound is larger in a fluid with higher temperature. However, for 
< A < the second term in the bracket in (1251) is positive which yields 
^‘s,QFT,i < ^%CFT < '^'s,QFT,r while Ti > Tr- Although this is a mathematically 
accepted solution, it is not physically relevant case. 

Moreover, plugging the expressions fl27|) into equations fllllllSp we will get the 
expression for the temperature of steady state as: 


T 


T. 



rA^ 


-i2(dH-l—A) 


+ * * • , 


(29) 


where Tc = is the temperature obtained for conformal case, and 


(1 - 

2d{d + l)(x + 1) 


[d{l + {d-6){l- , 


where we have dehned 



6 = 2 



(30) 


(31) 


Note that x > 1 since we have assumed Ti > T^. Examining the expression (122]) 
reveals that depending on the value of conformal dimension A the temperature T of 
the steady state starting with the initial perturbed fluids could be smaller or larger 
than the temperature Tc of the steady state starting with the initial conformal fluids. 
Indeed, 5 < 1 for A > {d+l)/2 which result in both terms in square bracket become 
negative. That is, in this regime r is negative and therefore T < T^. However, 5 > 1 
for A < (d + l)/2 which yields positive r and T > Tc. More precisely one gets 


T <Tr 


T>Tr 


for 

for 


d + 1 
2 

d — 1 


< A < d + 1, 
A d + 1 

< A < 


(32) 


2 2 

The above result is valid independent of the value of y, i.e., relative ratio of right 
and left temperature, however the temperature T deviates more from Tc when the 
difference between right and left temperature increases. 

Also one can obtain the boost velocity as: 


V = Vc. 


1 + 




T 


2(d+l-A) 


+ ■ • • , 


(33) 
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where Vc is given in equation fl2^ for the conformal fluid by replacing x ^ X 


r (d+l)^ x(x + l)(l + x 1 . 04 N 

[ 2d {x-l)[l-x-^)[x + d){x + d-^)\- ^ ^ 

The square bracket in the above expression is always positive then the sign of v is 
determined by the sign of r which we already discussed. Therefore from the expression 
([33]), one can deduce that the boost velocity is smaller or larger than the conformal 
case depending on the value of conformal dimension A, i.e., 


V > Vr_ 


V <Vr. 


for 

for 


d + 1 
2 

d-1 


< A < d T 1, 

A d + 1 

< A < 


2 2 

One can also obtain the shock speed in the left and right fluid as 


Ui — Ucj \ 1 -\- 


UiX^ 


2(d+l-A) 


T, 


Ur — "Uc, r I 1 T 




2{d+l-A) 


To 




+ ■■■ 


(35) 


(36) 


where Uc,i and Mc,r are given in equation fl2^ for conformal fluid by replacing x ^ X 
and 


v{d-l + 2x) (1 + d)V + a(l + d)x ^-a{d + x ^) + X ^-1 

1+d + dx-\x-m-^+x) 

^ z/(d - 1 + 2x~^) _ (1 + d)V + a{l + d)x^ - a(d + x) + X - 1 
1 + d (x-l)(d + x) 

It could be shown 



for 

d + 1 

<C A <C “h 1 , 


2 


for 

d-1 

A d + 1 


2 

2 ■ 


.(37) 

(38) 


(39) 


To check if the above solution is the unique solution of the problem, one needs to check 
the inequalities fll4p . From equation fl28p it is evident that for small perturbations, i.e., 
y/jnd+i-A ^ constraint in fll4p is always fulhlled. The second constraint 

is always satished for all Si and Sr for the conformal fluid and then it holds as well 
in our QFT model which is just a small perturbation around CFT. Therefore the 
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shock waves emanating at the contact point are stable and the results obtained for 
the steady state hold for all £i and £r- 

It might be interesting to compare the results of our QFT model in (1 + 1)- 
dimensions with that of [15]; note that the former is a perturbed CFT by an irrelevant 
operator while the latter is a CFT perturbed by a relevant operator. Let?s start with 
the left and right fluid which is described by fl27p with d = 1 and A = 1, i.e., 



(40) 


(41) 


The temperature and boost velocity of the steady state which develops in the inter¬ 
mediate region would be specified by and respectively, while the speed of 
shocks moving to the left and right are given by fl36|) . However, one would immedi¬ 
ately recover that the subleading terms in all the expansions vanish for d = 1 and 
A = 1 and obtain 


(42) 


T Tq , V Uc , Ul ^r,c 7 


for the steady state which is the same with that of conformal fluids. That is, even 
though the initial fluids that we start with on the left and right are perturbed CFTs 
given by EOS (jUj), the steady state developing in between left and right regions shows 
the same properties as one starts with conformal fluids as the initial states to start 
with. This is different from the results in [T5] for a steady state interpolating between 
left and right CFTs perturbed by a relevant operator where they observe deviation 
from the conformal case for the temperature, shock speeds and etc. 

We emphasis that expression flT2D is only true for steady state if the conformal 
dimension of the irrelevant operator is set to one in (1 -|- 1) dimensions. For any other 
allowed values of A one should use expansions 1321ES]) to express the properties of 
steady state and it definitely deviates from the conformal case. However, similar story 
happens in an arbitrary dimensions if one chooses an operator where A = (d -|- l)/2. 

3.3 Barotropic fluid 

Consider a barotropic fluid for which pressure is a function of energy density, i.e.. 


P = P{S ), 


(43) 
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where the energy density is only a fnnction of temperatnre, £ = £{T). We also 
assnme that the difference between left and right temperatnre is very small compared 
to either left or right temperatnre which implies 


£i — £r ^ £ ] 


Ur ? 


(44) 


with £i > £r. Now one can nse the Taylor expansion for pressnre around some To as 


P{£) — Pq + c^(£ — To) + ~ + ■ • ■ 


where we have dehned P{£q) = Pq and 


K = 


(PP 


(45) 


(46) 


is a characteristic parameter of the fluid and should satisfy the second inequality in 
ffTTjl . in order to have a unique shock solution to the Riemann problem with initial 
values £i for a; < 0 and T^ for x > 0. Now solving equations flllllldp perturbatively 
we will get the following expression for the energy density of the steady state fluid 

1 


T — To + 


16c? 




(47) 


where we have defined PIq = Pq + £q which is the enthalpy of the equilibrium state 
where T; = T^. Also, we have chosen To to be the mean energy of the left and right 
fluid, i.e.. 


Tn = 


£l + £r 


and 


AT = £] — £r 


(48) 

(49) 


Note that the first correction which appears in energy density is of second order. 
Indeed, this is not surprising as one would expect only even powers of AT appears in 
the energy expansion simply because the energy of the steady state should not change 
if we start by the initial condition where we have replaced T/ y-)- Sr- Furthermore, the 
boost velocity is given by 

n = ^AT + ---. (50) 

Here, we expect only odd powers of AT appears in the boost expansion since again 
the magnitude of the boost velocity should not change by replacing T/ y-)- T^, however 
one need to boost in the opposite direction. 
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The speed of sound in the left and right fluid to the hrst order in is given by 



'S 


(51) 


We did not include the second order term in the above expansion in order to avoid 
introducing a new characteristic parameter for the fluid. For k > 0 it is clear that 
Cs,l > Cs > Cs,r which is consistent with our assumption for the temperature that 
Ti > Tr. On the other hand for k < 0 we get Cs^r > Cg > Cg^i while we have assumed 
Ti > Tr- This is not what we usually expect in a normal physical system. Note that 
the speed of sound on left and right would be replaced under Si Sr- 

We can further obtain the speed of shock waves on the left and right fluid as 



(52) 


The bracket in the second term is always positive due to the stability constraint ffTT|) . 
therefore the left shock is subsonic compared to the speed of sound on the left while 
the right shock is supersonic with respect to the speed of sound on the right. 

If the characteristic parameter k > 0 then the second constraint in flTTl) is always 
satisfied and therefore the above solution is a unique solution. This is indeed a well- 
known fact, even for a classical fluid, that k > 0 is a sufficient condition for the 
existence and uniqueness of shock waves to the Riemann problem, known as Bethe- 
Weyl theorem However for an arbitrary equation of state with k < 0, it is 

still possible to have a unique shock solution for the Riemann problem as long as the 
inequality in ffTT|) fulhlled. 

4 Discussion 

We argued the possibility of performing the NESS when we bring two copies of sys¬ 
tems at different temperatures into a thermal contact in the framework of Bernard 
and Doyon [1], however, our systems enjoy more general EOS than the conformal 
fluid. The key feature in the BD set-up is that the steady state is not driven by 
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the temperature gradient. There are no external reservoirs, rather the initial sys¬ 
tems play the role of heat baths at sufficiently large scales. Furthermore, we used 
the holographic insight of [3] to describe the steady state as an intermediate Lorentz 
boosted state of the initial systems, interpolating between two asymptotic heat baths 
after long enough time. In this approach the formation of the steady state relies on 
the two single shock waves emanating at the point of contact moving in the opposite 
directions. Since both systems are assumed to be semi-inhnite then the claim is that 
there is no chance for the shocks to reflect back, forming a cascade and thermalizing 
the system. Although this seems to be sufficient condition to construct the steady 
state in the case of CFT, it is not enough for a system with a general EOS, e.g., 
barotropic fluid that we studied in this paper. In general, in order to perform the 
steady state, shock waves should not split otherwise thermalization may happen. This 
is due to the fact that if they split, they may move in the opposite directions and 
form composite waves leading to cascades and finally thermalizing the system. For 
a relativistic perfect fluid with general EOS the condition fll4p is sufficient in order 
to avoid splitting and to have stable shock solutions. Therefore the development of a 
steady state is guaranteed and one can obtain the properties of this NESS in terms 
of the characteristics of the fluid and the initial values, see equations flT7|l and fl50l) . 

Nonetheless, in our studies we only considered perfect fluids, for which viscosity 
is zero. It is interesting to see how the construction of the steady states would be 
affected in this set-up, if we move away from ideal hydro and take into account 
the viscosity. Since performing of the steady state depends on the propagation of 
shock waves as mentioned before, then we can turn around and ask how viscosity 
will affect the formation of shock waves at the point of contact. More precisely, one 
should investigate the possibility of shock solutions for the initial condition Riemann 
problem in viscous fluid. In fact, this question has been already studied in literature, 
e.g., see [IB] where the solutions to the relativistic Riemann problem for viscous fluid 
has been investigated numerically. As a result, by varying the ratio of shear viscosity 
to entropy density, i.e., 77 /s, from zero to infinity, a transition from ideal shock waves 
to viscous one has been shown; starting from ideal fluid with zero viscosity one obtains 
shock waves with zero width. By increasing rj/s the solutions with non-zero width will 
appear, the so-called viscous shocks. An upper limit for the tj/s has been estimated for 
which shocks can still be observed experimentally on the proper time scale. However, 
as one continues to increase the ratio r]/s above this upper limit, the free-streaming 
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will occur, i.e., shock solution is completely washed out. 
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